Separable sequences in Bianchi I loop quantum cosmology 
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In this paper, we discuss the properties of one-parameter sequences that arise when solving the 
Hamiltonian constraint in Bianchi I loop quantum cosmology using a separation of variables method. 
In particular, we focus on finding an expression for the sequence for all real values of the parameter, 
and discuss the pre-classicality of this function. We find that the behavior of these preclassical 
sequences imply time asymmetry on either side of the classical singularity in Bianchi I cosmology. 
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Loop quantum cosmology (LQC) [l| has been devel- 
oped as a symmetry reduced model of quantum geom- 
etry to test the viability of the latter as a theory 
of quantum gravity. Analogous to minisuperspace mod- 
els, LQC reduces the infinite degrees of freedom in the 
full theory down to a finite, quantum mechanical the- 
ory. This allows study of various models of cosmological 
interest, such as isotropic Friedman-Robertson- Walker 
(FRW) space-times. In particular, this has led to the 
emergence of a quantum gravity mechanism for infla- 
tion [2j, evolution through the classical singularity |j| 
or oscillatory universes [5j , as well as the development of 
effective semi-classical Hamiltonians |g. Developments 
in the toy model of LQC can also help direct inquiry in 
the full model 0. 

The quantum Hamiltonian constraint in LQC is a dif- 
ference equation, relating discrete eigenstates of the triad 
operators. The generic solution to a difference equation, 
however, has rather unphysical properties; in particu- 
lar, as the triad parameter is increased by discrete steps, 
the coefficients of the wave function expansion will al- 
ternate in sign as they increase in magnitude. This can 
be interpreted physically as quantum effects becoming 
important even in large-sized universes. At the present 
time, there is no physical inner product in LQC to elim- 
inate unphysical states with large quantum effects. This 
has led to the notion of pre-classicality ||, where states 
with oscillatory behavior far away from the singularity 
are eliminated from consideration. Generating function 
techniques 0, llf]| , useful in solving the multi-parameter 
recursion relations arising in the anisotropic models, can 
implement pre-classicality by placing boundary condi- 
tions on the function associated to the desired sequence. 

However, these methods have been used to find the 
solutions only for integer values of the parameters (al- 
though there exists a general solution in the isotropic 
case for all real values U see also ^3). This raises the 
question of the existence of solutions for non-integer val- 
ues in anisotropic models. Recent research has addressed 
this by treating each non-integer sequence distinctly |l2l | , 



rather than as part of a general solution for the recursion 
relation. In this work, we will find that obtaining the 
generating function allows us to find a single solution ap- 
plicable for all values of the parameters. We do this by 
separating the anisotropic recursion relations into prod- 
ucts of one-parameter equations. Models where this is 
possible include the Bianchi I system |9| and its locally 
rotationally symmetric (LRS) version [lOj; in addition, 
the Kantowski-Sachs model, relevant for the interior of 
the Schwarzschild black hole, is a similar case Be- 
cause the same one-parameter recursion relation appears 
in these contexts, it allows us to study the general prop- 
erties of wave functions as they pass from one side of the 
singularity to the other. Then we can examine their pre- 
classicality over the range of all triad eigenvalues, as well 
as differences in their behavior for positive and negative 
parameters. 

In particular, we look at the Bianchi I model, the 
anisotropic generalization of the spatially flat FRW cos- 
mology; a complete explanation of its quantization in the 
LQC context appears in 0|. As mentioned before, the 
quantized Hamiltonian constraint, acting on eigenstates 
of the triad operators, gives a difference equation relat- 
ing the coefficients s(mi , m,2, ms) of the wave function in 
the triad basis, ip = J] s(mi,m2, rns)\mi,m2, rn^). Note 
that the parameters rrik take real values, not just inte- 
gers. Because of the freedom to change the sign of two 
tetrad vectors, the coefficients s(mi, ra2, 77x3) in the full 
Bianchi I model must satisfy 

s(mi,TO 2 ,m 3 ) = s(—mi, —m 2 ,m 3 ) = s(—mi,m 2 , -m^) 
= s(mi, -m 2 , -m 3 ). 

There is a similar condition for the LRS case and its coef- 
ficients s(m, n), with m is the parameter in the direction 
of the rotational symmetry. This remaining gauge sym- 
metry allows us to let one of the parameters run from —00 
to 00, with the others restricted to non-negative values. 
Here, we choose the classical singularity for Bianchi I to 
be m3 = 0, and n — is used for Bianchi I LRS. 
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When the Hamiltonian constraint is quantized, 
we find that the coefficients s(mi,m 2 ,0) represent- 
ing the zero volume classical singularity fall out of 
the recursion relation. This is due to the dif- 
ference equation containing functions of the volume 
y(mi,ra 2 ,m 3 ) = {^j(p) 3 ^ 2 y/\mrm^m^\, where 7 is 
the Barbero-Immirzi parameter, and tp the Planck 
length. The actual difference equation is easier to deal 
with when we make the substitution i(mi, 7712, TO3) = 
y(2mi, 2m 2 , 2m3)s(2mi, 2m 2 , 2m 3 ), which gives rise to 
the conditions t(0, mi, 7713) = 0, etc. If we define the 
difference operator 8\ as 

5it(mi, 771.2,7713) = t{m\ + 1, to 2 , m 3 ) — t{m\ - 1, m 2 , m 3 ), 

and similarly for <5 2 and £3, the Hamiltonian constraint in 
the full Bianchi I model is given by the difference equa- 
tion 

[d(mi)5 2 ^3 + d(iri2)5iS 3 + d(m 3 )SiS2\t(mi, m 2 , m 3 ) = 0, 
i j (1) 

where d(n) 
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The factors of 



d(m k ) arise due to the various functions of the volume V 
existing in the constraint equation. In the LRS case, the 
Hamiltonian constraint becomes 

[d(n)S m 5 m + 2d 2 {m)5 m S n ]t(m,n) = 0, 

with 5 m and S n are the difference operators for the corre- 
sponding parameters, and again there are the conditions 
t(0,n) = t(m,0) = 0. When we write the equation in 
terms of a derived sequence i(m,n) = 8 m t{m,n) 1 then 
the relation for Bianchi I LRS becomes 



[d(n)S m + 2d 2 (m)S n ]t(m,n) = 0, 



(2) 



with the function d 2 (r7i) defined as 1/m when m ^ 0, and 
^2(0) = 0. Because of the switch from t(m, n) to t(m, n) 
in the LRS case, i(0, n) is no longer required to be zero. 

One method of solving the constraint relations and 
(0) is using a separation of variables technique. In partic- 
ular, if we look at the LRS case, we can assume the total 
sequence t(m,n) is the product a\(m)b\(n) of two one- 
dimensional sequences, with a separation parameter A. If 
we consider n, the parameter in the direction of the clas- 
sical singularity for Bianchi I LRS, as an "internal time" , 
then the sequence b\(n) describes the temporal evolution 
of the wave function. To solve the full recursion relation, 
the two sequences must satisfy 



5 m a\(m) = — a\(m), 



S n b\(n) 



-\d(n)b\(n), 



(3a) 
(3b) 



when m, n ^ (when m or n equal zero, the right hand 
sides of the equations are identically zero). The full 
Bianchi I case is similar - assuming the wave function 



coefficients t(m 1,1712,1x13) are the product of three one- 
parameter sequences c k (m k ), then each must solve 



8kC\ k {m k ) = \ k d(m k )c\ k (m k ), 



(4) 



with m k ^ and the separation parameters X k restricted 
by the relation AiA 2 + AiA 3 + A 2 A3 = 0. Again, choos- 
ing the singularity at m 3 = means that the sequence 
c\ 3 (7M3) gives the evolution of the wave function through 
the singularity. 

In the regime where n is small, the function d(n) ~ 
(l/2)<i 2 (n); far away from the singularity, the sequences 
solving the recursion relations 13a|) . (|3bl) and are es- 
sentially the same (up to a scaling of the separation con- 
stant A). This means that finding the solution for a par- 
ticular recursion relation is equivalent to finding one for 
all, at least for large values of the parameters. Closer to 
where m 3 or n are zero, the difference between the func- 
tions d(n) and c? 2 (m) appearing in the recursion relations 
will be more noticeable. Because the non-polynomial 
d(n) is more difficult to deal with, we shall first obtain 
the solution for the difference equation with rf 2 (m) , which 
has the same asymptotic values as one with d(n), then 
numerically evolve the sequence through the singularity 
using the function d(n). Another important point to note 
is the relative sign between the sequences. For example, a 
solution for a\(m) with a value A maps over to a solution 
to b\(m) with a separation constant of —A. Similarly, the 
condition AiA 2 + A1A3 + A 2 A 3 = for the full Bianchi I 
constraint means that at least one Afc must have a differ- 
ent sign from the others. 

We focus on the recursion relation (|3a|l for the sequence 
q,,\(t7?,) as a specific case, keeping in mind how to gener- 
alize this to the other recursion relations as discussed 
above. To solve this relation, we convert the problem 
of finding an appropriate sequence to that of finding the 
solution to a differential equation. This is a method that 
has been explored in detail elsewhere in the context of 
LQC so we only provide a summary here. One at- 
tempts to find a generating function F\(x), with a formal 
power series F\(x) = J2^ l=Q a\(m)x m , where the coeffi- 
cients a\(m) solve the recursion relation (|3a|l . Instead 
of working directly with the generating function F\(x) 
associated with the sequence, it turns out it is simpler to 
use the function G\(x), where F\(x) = a\(0) + xG\(x); 
this auxiliary function must satisfy 



d_ 

dx 



(l-x 2 )G x (x) 



2XG\(x) = a\(0), (5) 



with a\(0) — F\(Q) the value of the sequence at m = 0, 
and Ga(0) = a\(l). We will see that the values of a\(0) 
and A are related for pre-classical sequences. If A > 0, 
then any value of a\(0) will give prc-classicality; when 
A < 0, then sequences will be smooth on only one side 
of the singularity when a\(0) ^ and nowhere at all 
when a\(0) = (it is easily seen that the case A = 
corresponds to a constant sequence). 
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The situation a\(0) = enforces the boundary condi- 
tions at the singularity; the zero- volume part of the wave 
function drops out of the Hamiltonian constraint at this 
location. Solving the differential equation (JSJ for G\(x), 
and using this to find the generating function F\(x), we 
get 

F x (x) = G x (0)x(l + z) A - x (l - .t)-^ 1 . (6) 

The constant G\(0) acts as a scaling factor for the func- 
tion, and therefore the associated sequence. Thus, the 
solution we obtain for the recursion relation is essentially 
unique for a given value of A. Pre-classicality requires 
that there is no poles at x = —1; looking at the Taylor 



series coefficients of, e. g., (1 + x)^ 1 = 1 — x + x 2 + ■ ■ ■ 
shows that it has undesirable oscillations. This results 
in the condition A > 0. Once we have our generating 
function solution Jfyl, we can use the binomial theorem 
for each of the factors, and write out its series expansion. 
Then we read off the coefficients of the powers x m to get 
the solution to the recursion relation: 

^< m )^,(o)|Viy( m A _-/ 1 )(- A r 1 )^ 

This can be re-written in terms of hypergeometric func- 
tions as 



a\ >0 {m) = G\(0) 





i) 


\m — 





2 F 1 (l-m, A+l; A-m+1; -l)-(-l) 7 

I — 



A - 1 
-1 



-A-1 

m 



3F 2 (1, 1, A+m+Um+l, A+l; -1) 

m 



Even though this is coming originally from a Taylor 
series with positive integer values of m, we can use the 
formula (0 as the definition of the sequence for all val- 
ues of m, including the non-integers. An example of this 
exact solution is given in Figure Q for the value A = 3/2. 
The sequence smoothly passes through the classical sin- 
gularity Tfi — 0, and its behavior noticeably changes as it 
does so. Moving away from the singularity towards large 
positive m, the solution increases without limit; travel- 
ing from the singularity towards large negative m, the 
sequence undergoes oscillatory decay. Because these are 
the coefficients of the wave function describing the model, 
we can see that it is much more likely to find the Universe 
in the positive m regime. 

However, the recursion relation (|3a|l is only an ap- 
proximation to the usual case of when the singularity 
is crossed, such as the "temporal" sequences b\{n) for 
Bianchi I LRS, and ca 3 (777.3) for full Bianchi I. Exam- 
ining the behavior of these sequences as they cross the 
singularity will give insight into the generic properties 
of the complete sequences t(m, n) and t(mi, 7712,7713), re- 
spectively, for these types of models. Although we now 
are studying separable sequences where the prime appli- 
cation is to the sequences b\(n) and ca 3 (777,3), we will 
continue to use the notation a\ (m) as a generic name for 
all the solutions we find. 

As mentioned previously, the non-polynomial nature 
of the function d(n) appearing in the recursion relation 
is difficult to deal with using generating function tech- 
niques. We start with the exact values obtained by the 
definition 10 for large (non-integer) positive values of m 
and evolve numerically using equation Q with d(m); for 
the results shown in Figure [3 we have chosen the spe- 




FIG. 1: Graph of the exact solution a\(m) given in terms of 
hypergeometric functions by J7J for A = 3/2. For positive m, 
the sequence increases without limit. For negative m, there is 
an oscillatory decay in the wave function, representing a re- 
gion where there is a small probability of finding the Universe, 
compared to the positive side of the classical singularity. 



cific parameters A = 3/2 and in = k + 3/4, where k 
is any integer. Notice that, as the sequence passes the 
origin, the slope of the line changes, but there is no de- 
cay. Due to the invariance of the difference equation as 
777 — > —777, the sequence a(k + v), where k is an integer 
and v G [0, 1), is proportional to a(k+(l — is)), but the ra- 
tio of the two varies with v. Because of this dependence, 
the sequences parametrized by v should be thought of as 
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independent, rather than part of a unified solution for 
all to. The sequences have a deflection because of the 
behavior of d(m), differing greatly for \m\ < 1/2 from a 
simple 1/m dependence. This does not happen for the 
integer to sequences, since d(0) = 0, and the |m| < 1/2 
region is avoided. These results, and the sequence seen 
in Figure |2 are similar to those present in discrete quan- 
tum gravity 0]. Because we are seeing the same type 
of behavior for very different discretized version of loop 
quantum gravity, it is easy to suppose that these results 
are generic. Even though there may be modifications 
to the process of quantization - such as the use of a self- 
adjoint constraint in group averaging techniques - we can 
expect to see the same flavor of states as found here. 




FIG. 2: Graph of the numerically evolved solution a\(m) for 
A = 3/2 and using the recursion relation 13b^ with the func- 
tion d(m) and discrete parameters m = fc + 3/4, k an integer. 
Initial data is chosen for large positive m using the exact 
solution J7j, and evolved using the recursion relation in the 
negative direction. Note the large difference in slope between 
the two sides of the graph for both sequences. 

If we are focusing on one-dimensional sequences where 
the parameter is strictly non-integer, then we have not 
exhausted all the possibilities. We can consider the case 
when a\(0) ^ 0, which unfortunately makes the equa- 
tion JSJ for the generating function more complicated. 
There are two cases where sequences of this nature may 
arise. First, recall in the Bianchi I LRS model, when we 
simplify the Hamiltonian constraint under the definition 
t(m,n) = S m t(m,n), there is no restriction on the val- 
ues of t(0,n); however, because of the gauge conditions 
t(—m, n) — t(m, —n) there is no need to look at the so- 
lution for negative values of m. Another possibility is re- 
cent work |ll| using group averaging, where it was found 
that the physical Hilbert space consisted of sequences 
which were non-zero only at discrete values of the triad 
parameter to. Thus the sequence is typically zero at the 



classical singularity, yet if we find a pre-classical solution 
for all values, we can restrict this to only certain values 
and avoid the singularity in this manner. 

When A > 0, pre-classicality gives no conditions on 
the sequence, and the results are similar to the a\(0) = 
case - the solutions grow without bound for increasing 
values of the parameter m. However, when A < 0, there 
is a restriction relating the two values a\(0) and Ga(0); 
when this is used, the associated generating function is 
given by [lf1 | 

2 A x(l -x)-*- 1 



1- 



1-A 



2 F 1 (l-A,-A;2-A;(l+x)/2) 



Up to scaling by a\ (0) , the sequence is fixed by the value 
of A. As a simple case, we consider A = — 1, which gets rid 
of the (1 — x)~ x ~~ 1 factor in F\(x), with accompanying 
simplification in the Taylor series. Following the same 
procedure as with the a(0) = situation, we get 



a_i(m) 

w 



a-i 



= (-1)' 



(l-2mln2)-2 



J? (m 

3=0 



3 1) 



(8) 

This can be written as a complicated formula in terms of 
hypergeometric functions for all to, which we will not give 
here. However, this function is not defined at negative in- 
tegers, so we only use it for positive to, then numerically 
evolve for negative m with the recursion relation 13a(l . 
Because of the (— l) m factor, the function will be com- 
plex for non-integer to, but the ratio a_i(TO + l)/d— i(to) 
is real, so we can use the modulus of the function when 
finding the sequence values for positive to. The sequence 
gives behavior quite different from the a\(Q) = situ- 
ation as the singularity is passed. This can be see in 
Figure |3 where the oscillations on the left hand side are 
huge compared to the smooth evolution of the curve on 
the right hand side (too small to be seen in the graph). 
This result is generic for a\(0) ^ and A < 0; sequences 
that are pre-classical for positive values of to display large 
oscillations as the parameter becomes negative. The se- 
quences are unable to match the symmetry of the recur- 
sion relation, and therefore are only smooth on one side 
of the singularity. 

To summarize, we have examined separable solutions 
for the difference equations that appear as Hamiltonian 
constraints in LQC, and shown how to obtain them for 
all triad parameters, not just integer values. The large- 
scale behavior of these sequences is categorized by the 
separation constant A and value a a (0) at the singularity. 
In particular, for A < 0, the possible sequences displaying 
pre-classicality are very limited. When a\(0) ^ 0, the se- 
quences are smooth only on one side of the singularity; 
they do not exist at all for a\(0) = 0. On the other hand, 
for positive A, pre-classical sequences abound. There ex- 
ists a one-parameter family of solutions when a\(0) = 0, 
and for a\(0) ^ 0, the sequences are defined by their 
values of a\(0) and G\(0). 
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FIG. 3: Graph of the numerically evolved solution a\(m) for 
A = — 1, a(0) = 1 and using the recursion relation 13al with 
the function d,2{m) and discrete parameters m = fc + 3/4, k 
an integer. Initial data is chosen for large positive m using 
the exact solution ©, and evolved using the recursion rela- 
tion in the negative direction. Unlike the previous graphs, the 
sequence becomes oscillatory with greatly increasing magni- 
tude when m < 0. For positive m, the sequence is gradually 
decreasing, but much smaller magnitude (so it is not visible). 

These results have implications for the relationship be- 
tween pre-classicality and the large volume limit. As 
discussed using solutions of the Wheeler-de Witt equa- 
tion |lfj| . to assemble an arbitrary semi-classical profile 
for the wave function at large volumes, one needs a com- 
plete basis of functions made up of all possible values of 
A. For example, in Bianchi I LRS, when the temporal 
sequence b\(n) is restricted to integer sequences, then 
pre-classicality requires A < (because of the relative 
sign between the recursion relations). This gives spatial 
sequences a\(m) that decay to zero with increasing m. 
Non-integer sequences b\(n), which do not pass through 
the origin and hence are not fixed to zero there, can get 
around this limitation. This allows a total sequence made 
of sums across all A, so that any large volume behavior is 
feasible. The flip side is that this pre-classicality occurs 
only on one side of the singularity when A > (again, 
because of the sign flip for b\ (n)). 

Whether solutions to cosmological models should be 
pre-classical on both sides of the classical singularity or 
not is open for discussion. If one feels that in the distant 
"past" before the Big Bang, physics at scales compara- 
ble to the current scale of the Universe would necessar- 



ily be describable by partial differential equations, then 
pre-classicality on both sides of the singularity would be 
a physical requirement. In that case, the homogeneous 
models we have discussed here appear to lack the variety 
one would desire its pre-classical solutions to have. This 
would reinforce known results from the integer sector of 
the theory 0, [^}. On the other hand, if one insists on 
pre-classicality only on one side of the singularity, then 
there is no such issue. Both cases results in an inter- 
nal time asymmetry crossing the classical singularity - 
either the scale of the wave function changes (so physical 
observables such as volume vary with time at different 
rates), or the Universe descends into a quantum foam 
with no classical behavior at all. 
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